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1 Introduction 



Plasmodium is a vegetative stage of acellular slime mould Physarum polycephalum [40] 
The Plasmodium is a syncytium, a single cell with many nuclei, which feeds on mi- 
croscopic particles. During its foraging behaviour the Plasmodium spans scattered 
sources of nutrients with a network of protoplasmic tubes. The protoplasmic net- 
work is optimised to cover all sources of food and to provide a robust and speedy 
transportation of nutrients and metabolites in the Plasmodium body. The Plas- 
modium's foraging behaviour can be interpreted as computation: data are repre- 
sented by spatial configurations of attractants and repellents, and results of com- 
putation by structures of protoplasmic network formed by the Plasmodium on the 
data sets (33)35|6j . The problems solved by Plasmodium of P. polycephalum include 
shortest path [33135], implementatiton of storage modification machines |2j, Voronoi 
diagram [38j, Delaunay triangulation logical computing [ 44.4] . and process al- 
gebra [39]; see overview in [6]. 

Slime mould based computers are programmed using attractants and repellents. 
While reaction of P. polycephalum to repellents, e.g. illumination- [32], thermo- 
|43|29| and salt-based repellents [5], is quite straightforward the Plasmodium's be- 
haviour in presence of chemo-attractants still remains a challenging topic of re- 
search [I4] . 

Studies on attractants of P. polycephalum can be traced back to early 1900s, with 
interest usually reignited every 20-30 years, see overview of pre-1960s works in |9|. A 
common logic behind explaining attraction of Plasmodium to various substances was 
that slime mould prefers substances with potentially high nutritional value, partic- 
ularly those high on carbohydrates. Laboratory experiments in 1970s demonstrated 
that Plasmodium is attracted to glucose, maltose, mannose and galactose [9f25] . Re- 
cently there was a renewed interest in studying relations between chemo-attractants 
and their nutritional values [H] . Molecular mechanisms of carbohydrate recognition 
are proposed in [26j . 

Already in late 1970s it became clear that nutritional value is not the only pre- 
requisite for a substance to be a chemo-attractant for P. polycephalum [23|[M] . It 
was experimentally proved that the slime mould is attracted by peptones |12|9j , 
aminoacids phenylalanine, leucine, serine, asparagine, glycine, alanine, aspartate, 
glutamate, and threonine |11|23|30 ]. There are reports that Plasmodium is indif- 
ferent to sucrose, fructose and ribose [9 25j, and repelled by sucrose [46j and tryp- 
tophan [30j. Moreover, some of the chemo-attractants may inhibit Plasmodium's 
movement when the Plasmodium comes into direct contact with them, e.g. galactose 
and mannose attract Plasmodium [9,25j but also inhibit Plasmodium's motion [13] . 

While perfecting our concept and experimental implementations of Physarum ma- 
chines — programmable amorphous biological computers experimentally imple- 
mented in Plasmodium of P. polycephalum [6J we tried to find a universal way of 
programming plasmodium activity. A universal programming substance is the on 
that can be used to immobilise the plasmodium when necessary but also set the 
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Fig. 1. Photographs of a typical experimental setups, (a) Plasmodium propagates 
towards Kalms Sleep tablet, (b) Plasmodium envelopes a compact configuration of 
Kalms Sleep tablets. 



Plasmodium in a complex pattern of motion when required. We found that the 
herbal tablets with alleged sleep-inducing properties — Nytol, Kalms Sleep and 
Kalms Tablets — are universal programming substances. The paper presents out- 
comes of our studies and illustrates basic routines of programming the Plasmodium's 
behaviour. 



2 Laboratory methods 



We cultivated P. polycephalum in plastic containers, on paper towels sprinkled with 
still drinking water and fed with rolled oats. For experiments we used polystyrene 
Petri dishes (round, diameter 120 mm, and rectangular 120 x 120 mm), and 2% agar 
gel (Select agar, Sigma Aldrich) as a non-nutrient substrate and 2% Corn Meal Agar 
(Fluka Analytical, Sigma Aldrich) as a nutrient-rich substrate. 

We represented attracting domains with three types of sleep-inducing herbal tablets: 
Kalms Tablets and Kalms Sleep (G. R. Lane Health Products Ltd, Gloucester GL2 
OGR, UK) and Nytol (Brunei Healthcare Manufacturing Ltd, Derbyshire, DEll 
OBB, UK). We placed an oat flake colonised by Plasmodium of P. polycephalum at 
the fringe of a Petri dish and distributed few half-pills of Kalms Tablets, Kalms 
Sleep or Nytol in the central area of the dish (Fig. [T]). Images of Plasmodium are 
recorded by scanning Petri dishes in Epson Perfection 4490 (thus most images are 
views of Petri dishes from below). Photos are taken using FujiPix 6000 camera. 
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Fig. 2. Snapshots of Plasmodium propagating on a nutrient substrate with three 
Kalms Sleep (a) and Kalms Tablets (b) tablets. Snapshots are taken 20 h after 
inoculation of the Plasmodium in the centre of a Petri dish. 

3 Laboratory experimental results 

Finding 1 Herbal sleep-remedy tablets are treated by Plasmodium of P. polycephalum 
as obstacles when the Plasmodium grows on a nutrient- substrate. 

Being inoculated on a nutrient substrate Plasmodium of P. polycephalum prop- 
agates omni-directionally, in a fashion similar to an excitation wave in excitable 
medium [3J. When plasmodium approaches domains occupied by tablets it avoids 
close contact with agar gel nearby the tablets. Thus the Plasmodium's growth-front 
splits onto several wave- fragments (Fig. [2]). Depending on a physiological state of 
the Plasmodium and contamination of its substrate by other moulds and bacte- 
ria the plasmodial wave-fragments either propagate as isolated fragments or merge 
later in a single wave- front. The plasmodium reacts to the tablets as if they were 
sources of repellents, e.g. inorganic salts |46|5|6| . We did not observe any attraction 
of Plasmodium towards tablets in our experiments with nutrient agar gel. 

Finding 2 Herbal sleep-remedy tablets are treated by plasmodium of P. polycephalum 
as attractants when the plasmodium grows on a non-nutrient substrate. 

On a non-nutrient agar gel plasmodium propagates as localised excitation wave- 
fragments do in sub-excitable Belousov-Zhabotinsky medium [3J. It starts growing 
in one directions, usually follows gradients of attractants and repellents, humid- 
ity and illumination. Each growing 'cone', or active zone, the analog of travel- 
ling localisation, has a certain degree of autonomy yet, coordinates its strategy of 
space exploration with other parts of the plasmodium. In laboratory experiments we 
demonstrated that the plasmodium growing on a non-nutrient substrate is strongly 
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(c) (d) 

Fig. 3. Interaction of the Plasmodium growing on a non-nutrient substrate with a 
single Kalms Sleep tablet, (ab) Two stages of enveloping the tablet by Plasmodium, 
scanned Petri dishes, (c) Scheme of principal zones in plasmodium-tablet complex, 
(d) Photograph of the feeding zone, x40. 

attracted to the tablets but never comes into direct contact with the tablets or a 
substrate under the tablets. Further we are discussing experiments with the Plas- 
modium growing on a non-nutrient substrate. 

Let us consider experiments with a single tablet (Fig. |3]). Being inoculated at some 
distance from a tablet the Plasmodium propagates directly towards the tablet. On 
approaching the tablet, 2-4 mm away form the tablet, the Plasmodium generates 
two branches (Fig. |3^). One branch travels around the tablet clockwise another 
anti-clockwise. Eventually these two propagating fronts meet at the site opposite 
the site of the Plasmodium's initial branching, and merge (Fig. |3]3). A single tablet 
acts as a strong unequivocal attractive source in the Plasmodium's foraging space. 
After enveloping the tablet with its body the Plasmodium does not venture to any 
other domain of the substrate but stays virtually motionless. 
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(a) 12 h (b) 19 h (c) 29 h 




(d) 38 h (e) 47 h (f ) 

Fig. 4. Plasmodium, on a non-nutrient substrate, navigates around two Kalms Sleep 
(north part of Petri dish) and Kalms Tablet tablets (south part of Petri dish), 
(a-e) Snapshots of experimental Petri dish taken at various time intervals, indi- 
cated in hours from the moment of inoculation, (f) Scheme of the Plasmodium 
propagation: tablets are grey discs, inoculation site is the star and the trajectory of 
Plasmodium's propagation is shown by arrows. 

We can classify plasmodium-t ablet complex on the major transport artery, feeding 
zone and exploratory zone (Fig. [s]^). Major transport artery is a pronounced pro- 
toplasmic tube encircling the tablet. It separates Plasmodium's body on two parts: 
feeding zone and exploratory zone. Feeding zone (Fig. |3]l) consists of numerous 
branching processes staying at the edge of 'no-go' zone, proximal to the tablet's 
domain not occupied by Plasmodium. 

Finding 3 A single herbal tablet acts as a fixed attractor for Plasmodium of P. 
polycephalum while a sparse group of several herbal tablets may act as limit cycle. 

If two or more tablets are present in a Petri dish and 'no-go zones' of the tables 
do not intersect, the Plasmodium travels between the tablets in a periodic fashion. 
Here we consider three most common situations. Plasmodium travelling along a limit 
cycle around two tablets is shown in Fig. |4j Kalms Sleep tablet is positioned at the 
north part of Petri dish and Kalms Tablets tablet at the south part. Plasmodium 
is inoculated in the centre of the Petri dish. In first 12 h after inoculation the 
Plasmodium propagates north and encounters Kalms Sleep tablet (Fig. |4^). The 
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Plasmodium makes anti-clockwise turn around the tablet (Fig.|4]3) and travels south 
until it reaches Kalms Tablets tablet (Fig. The Plasmodium travels around 
the tablet anti-clockwise (Fig. [4]i) and heads north towards Kalms Sleep tablet 
(Fig. [4^). If there was no accumulation of Plasmodium's metabolites in agar gel the 
Plasmodium would travel around the tablets non-stop (Fig. [4]F). 



In example shown in Fig. |4] plasmodium keeps its structural integrity during cir- 
cumnavigation. Active zone remains connected to the rest of Plasmodium's body 
with a network of protoplasmic tubes. Let us now consider experiment where Plas- 
modium splits into several independent plasmodia, each one with its own active 
zone (Fig. [5]). While describing Plasmodia's behaviour we will refer to Nytol tablets 
numbered from 1 to 4 as in scheme Fig. [5^h. At the beginning of experiment an oat 
flake colonised by plasmodium is placed in north-east corner of a rectangular Petri 
dish. In 13 h after inoculation an active zone is formed. The active zone propagates 
towards tablet 1 and splits into two active zones (Fig. [5^). One active zone does not 
survive, another travels around tablet 2 and heads towards tablet 4 (Fig. [5]3). The 
active zone reaches tablet 4 and starts moving clockwise around the tablet by 52nd h 
of experiment (Fig.[5]3). After passing tablet 4 along its western side the plasmodium 
reaches tablet 3 and turns around the table 3 clockwise (Fig. [sjl). At roughly the 
same time the plasmodium still residing on the colonised oat flake emits new active 
zone which travels north of tablets 1 and 2 towards tablet 3 (Fig. [5]l). The active 
zone moving around tablet 3 and heading east-south-east collides with active zone 
travelling west around tablet 2. In the result of such collision two active zones merge 
and form a single active zone travelling south (Fig. [5^). This active zone reaches 
tablet 4, turns around the table clockwise and travels north (Fig.[5]F). The scheme of 
propagation is shown in Fig. [6^. In an ideal situation, the system plasmodia-Nytol 
tablets would form a carousel of active zones rotating around tablets 3 and 4 and 
supported by streams of plasmodia periodically emitted from the site of inoculation 
(Fig.|6}3). 



Example of complex navigation of plasmodia around conflguration of Nytol tablets 
is demonstrated in Fig.[7|. An oat flake colonised by plasmodium is placed in north- 
east corner of rectangular Petri dish. The plasmodium propagates towards tablet 
1 (see labelling of tablets in Fig. [7^) and splits into two active zones A and B 
(Fig. [7^b). In 40-50 h after inoculation these active zones become disconnected 
and travel as independent plasmodia (Fig. [7j3c). Plasmodium A travels north of 
tablet 2 and plasmodium B heads towards tablet 4 (Fig. At the same time 
additional active zone C is emitted from the oat flake colonised by plasmodium 
and propagates towards tablet 1 (Fig. [tJI). By 60th hour from the beginning of 
experiment plasmodium B makes almost complete clockwise rotation around tablet 
4 while Plasmodium A makes anti-clockwise rotation around tablet 2 (Fig. [tJI). 
Eventually plasmodia A and C collide, merge and propagate further as a single 
Plasmodium D (Fig. [7^). This plasmodium D rotates clockwise around tablet 3 
and Plasmodium B continues its clockwise rotation around tablet 4 (Fig. [7]F). See 
trajectories of active zones and plasmodia in Fig. [S) 
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(d) 68 h 




(e) 90 h 



(f) 113 h 



Fig. 5. Formation of the 'limit cycle' supported by streams of Plasmodium 
wave-fronts, on a non-nutrient substrate, (a-f) Snapshots of plasmodia propagating 
in the configuration of Nytol tablets. 
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(a) 



(b) 



Fig. 6. Schematic formation of the 'hmit cycle' supported by streams of Plasmodium 
wave-fronts, on a non-nutrient substrate, (a) Trajectories of plasmodia propagation, 
solid line indicates plasmodia propagation during first 60 h after inoculation, dotted 
line shows trajectory of Plasmodium emerged and developed during period 68-90 h, 
and double line reflects propagation during 90-113 h. (b) An idealised situation when 
active zones rotate around the tablets fed by streams of plasmodia. The schemes 
are derived from experimental data presented in Fig. [5j 

4 Numerical simulation 



A profile of Plasmodium's active zone growing on a non-nutrient substrate is iso- 
morphic to shapes of wave- fragments in sub-excitable media When active zone 
of P. polycephalum spreads, two processes occur simultaneously — advancing of the 
wave-shaped tip of the pseudopodium and formation of the trail of protoplasmic 
tubes. We simulate the chemo-tactic travelling of Plasmodium using two-variable 
Oregonator equations jl7|45j : 




The variable u is abstracted as a local density of Plasmodium's protoplasm and 
V reflects local concentration of metabolites and nutrients. We integrate the sys- 
tem using Euler method with five-node Laplace operator, time step At = 5 • 10~^ 
and grid point spacing Ax = 0.25, with the following parameters: = 00 — ?7/2, 
A = 0.0011109, 00 = 0.0766 for imitating plasmodium growth on a non- nutrient 
substrate and 00 = 0.066 for nutrient-substrate, e = 0.03, / = 1.4, q = 0.022. 

Parameters q and / are inherited from the Oregonator model of Belousov-Zhabotinsky 
medium [17 45j, is proportional to local concentration of attractants and repel- 
lents. The parameter 77 corresponds to a gradient of chemo-attractants emitted by 
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(a) 12 h 



(b) 37 h 





(c) 49 h 



(d) 61 h 





(e) 85 h 



(f) 109 h 



Fig. 7. Complex navigation of plasmodia in the configuration of four Nytol tablets, 
on a non-nutrient substrate, (a-f) Snapshots of experimental Petri dish. 
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Fig. 8. Scheme of complex navigation of plasmodia in the configuration of four Nytol 
tablets, on a non-nutrient substrate, derived from experimental data in Fig. [7| 

data planar points. Let P be a set of attraction sites P and x be a site of a sim- 
ulated medium then r]x = 2 • 10~^ — mmp^p{d{x^p) : ^{p) = True} • where 
3.1 • 10^ < 6 < 4.9 • 10^ and d{x,p) (for the simulated medium 400 x 400 sites) is an 
Euclidean distance between sites x and p. 



We imitate narrow feeding zones around tablets (Fig. [3pib) as follows. Let x be a 
centre of a tablet radius r > 0, we assume (j)o = 0.065 in any site y such that 
^ ^ d{y, x) < r ^ where is a the feeding zone's width. 



The medium is perturbed by an initial excitation, where a 11 x 11 sites are assigned 
1^ = 1.0 each. The perturbation generates a propagation wave- fragment travelling 
along gradient 77. Repellents emitted by sites of P show very limited diffusion. 
Therefore repellents can be regarded as impassable obstacles. 

To imitate formation of the protoplasmic tubes we store values of u in matrix 
L, which is processed at the end of simulation. For any site x and time step t if 
Ux > 0.1 and Lx = then Lx = 1. The matrix L represents time lapse superposition 
of propagating wave-fronts. The simulation is considered completed when propagat- 
ing pattern envelops P and halts any further motion. At the end of simulation we 
repeatedly apply the erosion operation [6], which represents a stretch-activation 
effect [22J necessary for formation of Plasmodium tubes, to L. The resultant proto- 
plasmic network provides a good phenomenological match for networks recorded in 
laboratory experiments. 

Our experimental finding on behaviour of Plasmodium on a nutrient substrate 
(Fig. [2]) is fully verified in Oregonator computer model. The tablets are treated 
as obstacles and no attraction of Plasmodium wave-fronts to the tablets is observed 
in computer simulations (Fig. l9]). 



Simulation of Plasmodium enveloping a single table is shown in Fig. 10 We generate 
a source of attractants at the south-west corner of a simulation space and perturb 
the medium in a 15 x 10 sites domain in the north-east corner of the simulation 
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(a) (b) 

Fig. 9. Simulation of Plasmodium propagation on a nutrient substrate with four 
tablets, (a) Initial configuration, tablets are shown as blue/dark grey washers and 
initial position of Plasmodium is a red/light grey ellipse, (b) Time-lapse snapshots 
of propagating Plasmodium's front. 



space (Fig. [lOpi). Active zones, visible as fronts of wave- fragments in Fig. 10 3-f 
travel along attractant gradients towards the tablet. On approaching the tablet 
the Plasmodium's active zones encounter repellents and therefore do not propagate 
close to the tablets but remain in the feeding zone (Fig. [lOj^de) and encapsulate the 
repellent domain (Fig. [To]F). When active zone growing clockwise collide with active 
zone growing anti-clockwise the zones merge (Fig. ^Of). 



Network of major protoplasmic tubes extracted from lapse-time images (Fig. [T0| ) is 
shown in Fig. [TT^ . The network satisfactory approximates the protoplasmic network 
recorded in laboratory experiments (Fig. [S^b). Summary scheme of Plasmodium's 
propagation towards and enveloping of a simulated tablet is shown in Fig. [TT) 3. 

Simulated periodic motion of Plasmodium's active zone around tablets is shown 
in Fig. [T2j The simulation results are in agreement with experimentally observed 
motion patterns (Figs. [4| [5]and[7|. We arrange two simulated tablets as shown in 
Fig. [12^, and offset initial perturbation domain (analog of Plasmodium's inoculation 
site) some distance away from a line passing through the tablets' centres. The tablets 
are arranged such that there is a space between the tablets exceeding double width 
of the feeding zones. 

When Plasmodium propagates close to tablet 1 the feeding zone splits in two parts, 
one part travels anti-clockwise another clockwise. Due to initial offset the zone trav- 
elling anti-clockwise annihilates and only active zone travelling clockwise continues 
its journey (Fig. [T2|3). When the active zone, travelling around tablet 1, approaches 
tablet 2 (Fig.[T2|3), it branches to a second active zone, which starts travelling clock- 
wise around tablet 2 (Fig. [T2|d). The active zone travelling clockwise around tablet 
1 then finds itself in a close proximity of both tablets. Therefore the active zone 
splits into two propagating patterns: one active zone continues travelling clockwise 
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inoculation 




tablet 



(a) 



(b) 




(c) 



(d) 



(e) 



(f) 



Fig. 10. Simulation of Plasmodium enveloping a single tablet on a non-nutrient 
substrate, (a) Initial setup, (b-f) Time-lapse images of propagating active zones. 
Zone of repellents is shown by blue/grey in (a). 
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Fig. 11. Simulation of Plasmodium enveloping a single tablet on a non- nutri- 
ent substrate, (a) Major protoplasmic network extracted from time-lapse images, 
(b) Scheme of interaction of active growth zones. Inoculation site is shown by star 
in (b). Trajectories of Plasmodium active zones' propagations are shown by arrows. 

around tablet 1 while second travels anti-clockwise around tablet 2 (Fig. [12^). 

Therefore at some stage we observe three active zones (wave-fragments) circum- 
navigating the tablets: one active zone running clockwise around tablet 1, and two 



zones running anti- and clockwise around tablet 2 (Fig. 12 3f). The active zones 
travelling around tablet 2 collide head-on (Fig. [T2^f) and merge. The active zone 
travelling around tablet 1 continues its travel and the cycle of splitting/branching 
and annihilation is repeated from the step shown in Fig. [l2b. 



The static image of active zones' trajectories (Fig. |13p i) and the reconstructed net- 
work of protoplasmic tubes (Fig. fl3b) are similar to Plasmodium shapes observed 



in laboratory experiments. See also a scheme of periodic movement in Fig. [T3p . In 
laboratory experiments periodic movements never last more than 3-5 repetitions. 
Usually at some stage resources are depleted, or a substrate is contaminated with 
metabolites, and the Plasmodium migrates to other parts of experimental arena. 



5 Possible applications of herbal tablets in computing with slime 
moulds 



In laboratory experiments and computer simulation we demonstrated that herbal 
tables may be in an untraditional ways of controlling spatial development of slime 
mould Physarum polycephalum. The slime mould is attracted to a single tablet, 
envelops a group of tightly packed tablets and circles between loosely positioned 
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inoculatfon 




(e) (f) 

Fig. 12. Simulation of the Plasmodium orbiting around two tablets on a non-nutrient 
substrate, (a) Initial setup, (b)-(f) Snapshots of wave-fragments, corresponding to 
tips of active zones, (g) Time-lapse image of propagating active zones. 
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(c) 

Fig. 13. Simulation of the Plasmodium orbiting around two tablets on a non-nutrient 
substrate, (a) Time-lapse image of propagating active zones, (b) Major protoplasmic 
network extracted from time-lapse images, (c) Scheme of periodic propagation of 
active zones. 

tablets. Possible active substances which could be responsible for such an unusual 
interaction between slime mould and herbal tablets are analysed in [7J. The herbal 
tablets is a cheap, user- friendly, and efficient alternative to existing techniques of 
controlling slime moulds with illumination- [32j, thermo- [43 _ 29j . and salt-based 
repellents [5], and carbohydrate-based attractants j9|25|14j . In present section we 
discuss two possible applications of sedative tables in slime mould based computing: 
approximation of planar shapes and implementation of logical gates. 



5. 1 Planar shapes 

A convex hull of a finite set P of planar points is the smallest convex polygon that 
contains all points of P (Fig. [14^). a-hull of P is an intersection of the complement 
of all closed discs of radius 1/a that includes no points of P [T5fT6] . a-shape is 
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(a) (b) 

Fig. 14. Planar hulls, (a) Example of convex hull, (b) Example of concave hull. 
Points of set P are empty discs. 



a convex hull when a ^ oo. With decrease of a the shapes may shrink, develop 
holes and become disconnected, the shapes collapse to P when a ^ 0. A concave 
hull is non-convex polygon representing area occupied by P. A concave hull is a 
connected a-shape without holes (Fig. [l4|3). We want slime mould to solve the 
following problem. Given planar set P represented by physical objects Plasmodium 
of P. polycephalum must represent concave hull of P by its largest protoplasmic 
tube. 

The fist algorithm of convex hull construction [19] was based on cognitive tactic 
techniques we use in our everyday's life. We select a starting point which is ex- 
tremal point of P. We pull a rope (anti-) clockwise to other extremal point. We 
continue until the set P is wrapped completely. The computation stops when we 
reach the starting point. Let we represent data points P by sources of attract ants 
only, e.g. by oat flakes (Fig. [15^). We place a piece of Plasmodium at some distance 
away from the set of points (Fig. [l5^). The Plasmodium propagates towards set 
P, colonises oat flakes (Fig. [Ts^b) and spans them with a network of protoplasmic 
tubes (Fig. fi5\p). No hull is constructed. When all data points P are colonised and 
spanned by protoplasmic network the Plasmodium ventures to explore the space 
around P and thus loses all chances of constructing any shapes of P. 

Proposition 1 P. polycephalum does not compute concave or convex hull of a set 
represented by attracting sources. 

Using repellents only, e.g. as a set of obstacles between attractants and inoculation 
site, will do no good because Plasmodium will just pass around the repellents and 
leave (see details in [6J, controlling Physarum with salt). The only solution would 
be to employ attractants to 'pull' Plasmodium towards planar set P and to use 
repellents to prevent Plasmodium from spanning the points of P (Fig. [l6]). Strength 
of repellents should be proportional to a and thus will determine exact shape of the 
constructed hull. This corresponds to original definition fT5] that a-hull of P s the 
intersection of all closed discs with radius l/a that contain all the points of P. 
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(c) 

Fig. 15. Slime mould spans oat flakes, (ab) Spanning oat flakes, representing P, 
by a network of protoplasmic tubes: scanned image (a) and its binarisation (b). 
(c) Scheme of the Plasmodium propagation, circles are points, or oat flakes, of P 
and star marks site of inoculation, arrows are protoplasmic tubes. 

Proposition 2 Plasmodium of P. polycephalum approximates connected a-hull with- 
out holes of a finite planar set, which points are represented by sources of long- 
distance attractants and short- distance repellents. 

When presented with a half-pill of the Kalms Tablets/Sleep the Plasmodium prop- 
agates towards the pill and forms, with its protoplasmic tubes, a circular enclo- 
sure around the pill. Such a unique behaviour of Plasmodium in presence of Kalms 
Tablets/Sleep indicates that a Plasmodium could implement Jarvis's Gift Wrapping 
algorithm [19], adapted to concave hulls, if points of P are represented by the pills. 

We tested feasibility of the idea in 25 experiments. All experiments were successful. 
In each experiment we arranged 4-8 half-pills (they represent given planar set P) in 
a random fashion near centre of a Petri dish and inoculated an oat flake colonised 
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\ 



Fig. 16. Proposed distribution of attracting and repelling gradients which may force 
Plasmodium to approximate a concave hull (arrows aiming towards set of discs are 
attractive forces, arrows originating in data points, discs, are repelling forces). 



by Plasmodium 2-4 cm away from the set P. A typical experiment is illustrated in 
Fig. [T7| In 12 h after inoculation plasmodium propagates towards set P and starts 



enveloping the set with its body and network of protoplasmic tubes (Fig. 17). The 



Plasmodium completes approximation of a shape by entirely enveloping P in next 



12 h (Fig. 17). The plasmodium does not propagate inside configuration of pills 



(Fig. [T7|). Scheme of the concave hull construction is shown in Fig. 18 



Configuration of P in Fig. 17a.-d favours approximation of a convex hull. If spatial 



configuration of points curves inwards then concave hull is approximated (Fig. 19). 



Typical hulls approximated by plasmodium of P. polycephalum in our experiments 
are shown in Fig. [20) Hulls constructed by plasmodium match their counterparts 
calculated by classical algorithms. 



Let us verify experimental laboratory results with computer simulation. To imitate 
formation of the protoplasmic tubes we store values of u matrix L, which is processed 
at the end of simulation. For any site x and time step t if Ux > 0.1 and Lx = 
then Lx = I. Gradients are chemo-attractants used in simulation are visualised in 
Fig. [2T^ . The matrix L represents time lapse superposition of propagating wave- 
fronts (Fig.[2T|3). The simulation is considered completed when propagating pattern 
envelops P and halts any further motion. At the end of simulation we repeatedly 
apply the erosion operation [6J (which represents a stretch-activation effect |22] 
necessary for formation of plasmodium tubes) to L. The resultant protoplasmic 
network (Fig. provides a good phenomenological match for networks recorded 
in laboratory experiments. 
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(c) t = 24 h (d) t = 24 h 



Fig. 17. Computation of concave hulls of point set P. Photographs (ac) and their 
binarised images (bd) taken 12 h (ab) and 24 h (cd) hours after Plasmodium's 
inoculation. 

5.2 Logical gates 

It was already demonstrated that Plasmodium can approximate Voronoi diagram 
and Delaunay triangulation, compute proximity graphs, approximate shortest path [6J 
Experimental implementation of concave hull computing by P. polycephalum [8\ put 
a logical closure in our designs of real-life prototypes of Physarum machines for com- 
putational geometry. Collision-based computing P could be yet another domain of 
unconventional computing where the unique interaction of slime mould with herbal 
tablets can be explored. Previously we have demonstrated [4J that due to inertial 
component of its movement active zones of P. polycephalum implement basic logi- 
cal gates in a geometrically constrained substrate. They propagate along channels 
cut from agar plate, one Plasmodium can block and deflect propagation of another 
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Fig. 18. Computation of concave hulls of point set P. Scheme of the Plasmodium 
interaction with P: points of P are circles; diffusing repellents, visible as black halo's 
surrounding pills in Fig. [iT^c, are shown by gray gradients; feeding boundary of 
protoplasm is marked by dotted line; major protoplasmic tube, which approximates 
Ci7(P), is shown by solid line. 




Fig. 19. Binarised image of Plasmodium approximating concave hull of P. 

Plasmodium PJ. Mechanics of interaction of Plasmodium with sedative tablets of- 
fer an opportunity to implement plasmodium-based logical gates in a free space, 
architect urelss, medium. 

Proposition 3 Using herbal tablets we can implement functionally complete set of 
logical functions: conduction, negation and disjunction. 

In a framework of collision-based computing [Ij, propagating active zones of the 
Plasmodium, e.g. zones A-D in Fig. [7| are considered as 'quanta' of information: 
a presence of a zone in a given space domain instantiates Boolean Truth and 
absence of the zone in the domain Boolean False. When two active zones collide 
they change their velocity vectors, e.g. deflect or merge. Post-collision trajectories 
of active zones represent Boolean values of logical operation computed in the zones' 



21 



Adamatzky A. Mathematical and Computer Modelling 55 (2012) 884-900. 




Fig. 20. Contours of convex and concave hulls constructed on small sets of planar 
points are shown. The contours are extracted in laboratory experiments with P. 
polycephalum. 



collision. 

Consider scenario shown in Fig. [7]le. Active zone A collides with active zone C. 
They fuse in new active zone D. If zone A (C) was absent then zone C (A) would 
continue along its original trajectory (Fig. [22^). Assuming zone A represents logical 
value of a Boolean variable x and zone C value of a variable y the colliding zones 
implement the gate (x^y) (xy^xy^xy). Trajectory of propagation of undisturbed 
(i.e. not colliding with zone A) zone C represents xy, undisturbed trajectory of zone 
A represents xy. Zone D, produced in the result of the collision between A and C, 
represents conjunction xy of input variables x and y. 

If we assume that input ^ is a constant Truth then negation gate will be as in 
Fig. [22) 3. This is a two-output gate, one output is a logical negation and other 
output is identity. A construct of a disjunction gate is shown in Fig. [22| 3. A func- 
tionally complete logical system is universal. Thus we provided a sketchy proof that 
a plasmodium-tablets system is a universal computing device. 

All these results are rather indicative. More experiments — in laboratory and com- 
puter simulation — will be required to develop a 'universal programming pill', which 
will combine all necessary components to program advanced computation with acel- 
lular slime mould P. polycephalum 
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(c) 

Fig. 21. Simulating plasmodium approximation of convex hull of planar set P in 
Oregonator model, (a) Gradients of repellents generated by points of P. (b) Time 
lapsed images of propagating plasmodium patterns, (c) Structure of protoplasmic 
tubes developed. 

6 Discussion 



We developed computer model imitating behaviour of slime mould P. polycephalum 
in presence of herbal sedative tablets. In laboratory experiments and numerical 
simulation we demonstrated that the herbal tablets Nytol, Kalms Sleep and Kalms 
Tablets play a role of fixed attractors and limit cycle in the space-time dynamics 
of the Plasmodium of P. polycephalum. The plasmodium reacts to presence of the 
tablets in its environment in a particular manner, drastically different from how 
the Plasmodium reacts to conventional sources of nutrients, esp. carbohydrates. 
When presented with a configuration of rolled oats, or oat flakes, the plasmodium 
propagates towards the flakes and spans them with a network of protoplasmic tubes 
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(b) 



x+y 

(c) 



Fig. 22. Scheme of collision-based gates implement able in plasmodium-tablets sys- 
tem, (a) (x, y) {xy, xy, xy) gate can be implemented in interaction between active 
zones A and C as illustrated in Fig. [7]le, solid arrows are incoming trajectories of 
active zone A (variable x) and active zone C (variable dotted arrow is outgoing 
trajectory of active zone dashed arrows represent undisturbed trajectories of 
zones A and C . (b) Negation, (c) Possible implementation of conjunction. 




(a) 



(b) 



Fig. 23. (a) Plasmodium spans oat flakes by a network of protoplasmic tubes, 
(b) Plasmodium occupies oat flakes (shown by arrows) but does not come into 
direct contact with Kalms Tablets, In all experiments illustrated Plasmodium grows 
on a non-nutrient substrate. 



(Fig. |23p.). In a mixed configuration of oat flakes and herbal tablets, oat flakes are 
still used as sources of nutrients however a presence of the oat flakes does not change 
overall pattern of Plasmodium's foraging behaviour (Fig. [23p). 



When given a choice between a classical source of carbohydrates, oat flakes or 
pasta/bread or even honey, and the herbal tablets the Plasmodium prefers herbal 



tablets in over 70% of experiments [7j, see example in Fig. 24 a.. The Plasmodium de- 
tects location of herbal tablets even when they are isolated from agar gel substrate. 
Thus in experiments illustrated in (Fig. [24)3c) we placed tablets on polyethylene 
pads, which prevented diffusion of substances from tablets to agar gel. Thus plas- 
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(a) 



(b) 




(c) 



Fig. 24. (a) Plasmodium is attracted to Kalms Tablets tablet placed on plastic pad 
(right) and does not propagate towards oat flake (left), (b) Plasmodium detects 
Kalms Tablets via airborne gradients of chemo-attractants and attempts to ap- 
proach two tablets, (c) scheme of the experiment (b). In all experiments illustrated 
Plasmodium grows on a non-nutrient substrate. 

modium was unable to rely on gradients of chemoattractants in its substrate but only 
the chemo-attractants in the air. In all such experiments Plasmodium successfully 
detects the tablets and propagates towards them. As we see in example Fig. [24|3C 
the Plasmodium detects location of tablet 1, propagates towards this tablet and 
tries to climb the polyethylene pad to get closer to the tablet. The Plasmodium fails 
and then it propagates towards tablet 2 and makes few more attempts to mount 
the tablet's plastic pad. 



Based on our experiments we can assume that the herbal tablets contain at least 
two types of chemical substances: one substance plays a role of long-distance chemo- 
attractants, another substance is a short-distance repellent. Amongst six herbal 
extracts present in the tablets [27j, only Valerian roots (Valeriana officinalis) and 
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Hops {Humulus lupulus) extracts are present in all three types of tablets. Thus we 
can propose either Valerian or Hops is a long-distant attractant for P. polycephalum. 
To check which one, we undertook an extensive range of experiments with dried 
plants and found that with probability 0.73 Plasmodium of P. polycephalum prefers 
dried Valerian root to dried Hops leaves, see details in [7J. 

What exact substance in Valerian attracts Plasmodium remains unclear. Valerian 
contains hundreds of identified components including alkaloids [42j, volatile oils [18], 
valerinol [20], and actinidine [21j. Isovaleric acid and actinidine are identified in 
anal gland secretion of Iridomyrmex nitidiceps ant, and isovaleric acid is considered 
to be a distress indicator [lOj. Possibly there is a link between the said findings, 
and isovaleric acid combined with actinidine could be considered as pheromones of 
P. polycephalum (we could cite only relevant studies conducted on cellular slime 
moulds [28 l37|31j ). By allowing the Plasmodium to choose between Valerian and 
Catnip (Nepeta cataria) we ruled out actinidine, because in all experiments P. 
polycephalum prefers Valerian to Catnip. 

Amongst components of the tablets' filling and coating only stearic acide and mag- 
nesium stearate are reportedly present in all three tablets Also Kalms Sleep 
and Kalms Tablets tablets contain substantial (up to 160-200 mg [27]) amount of 
sucrose. Influence of stearic acid and magnesium stearate on slime mould is not 
studied, and we were unable to find any published data. Sucrose would be main 
candidate on a role of short-distance repellents, responsible for the formation of a 
'no-go zone' (Fig. [sj^) and, possibly, formation of feeding zone in the region with 
low concentration of a sucrose. Reports about influence of sugars are on the Plas- 
modium are somewhat contradictory. Some studies suggest that Plasmodium of P. 
polycephalum is indifferent to sucrose, fructose and ribose [9]25j while others demon- 
strate that Plasmodium is actually repelled by sucrose [46j . Galactose and mannose 
reportedly inhibit Plasmodium's motion, in situations of direct contact, yet they do 
not act as long-distance repellents [13j. 

Despite the fact that exact mechanism of the herbal tablets influence on the Plas- 
modium of P. polycephalum remains unclear, our present findings convincingly 
show that the herbal tablets containing Valerian extract (and possibly Hops) could 
be efficiently used as control stimuli in programming and controlling behaviour 
of Physarum machines [6J. The herbal tablets is a cheap, user- friendly, and effi- 
cient alternative to existing techniques of controlling slime moulds: illumination- 
[32j, thermo- |43|29 j. and salt-based repellents [5j, and carbohydrate-based attrac- 
tants [ 9l25ll4j . 

What are potential application domains of our results? Indeed several important 
questions on sensitivity of Plasmodium to active substances present in sedative 
tablets must be answered in a way of rigorous biological experiments. Exact mecha- 
nisms of slime mould's chemo-attraction are still to be discovered and Plasmodium 
preferences to various chemical components are to be classified. In the framework of 
biological computing our results provide a viable alternative to existing approaches 
and techniques of programming experimental laboratory biological substrates. A 
single table can cause a wide range of behavioural patterns in slime mould and 
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therefore — ideally — the whole program of slime mould based computer can be 
encoded in a single tablet. Slime mould of P. polycephalum could also play a role 
of an experimental laboratory model in design of amorphous soft-bodied robots 
and their application of substance transportation in confined environment and drug 
delivery inside human bodies. Our findings on control using combinations of long- 
distance attracts and short-distance repellents will allow for a precise navigation 
and control of such amorphous robots. 
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